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ABSTRACT 

We use a wide range of observations to constrain cosmological models possessing a 
significant asymmetry in the lepton sector, which offer perhaps the best chance of 
reconciling a critical-density Universe with current observations. The simplest case, 
with massless neutrinos, fails to fit many experimental data and does not lead to 
an acceptable model. If the neutrinos have mass of order one electron-volt (which 
is favoured by some neutrino observations), then models can be implemented which 
prove a good fit to microwave anisotropics and large-scale structure data. However, 
taking into account the latest microwave anisotropy results, especially those from 
Boomerang, we show that the model can no longer accommodate the observed baryon 
fraction in clusters. Together with the observed acceleration of the present Universe, 
this puts considerable pressure on such critical-density models. 

Key words: cosmology: theory - large-scale structure of the Universe 
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1 INTRODUCTION 

The recent use of the magnitude-redshift relation of type la 
supernovae to infer that the present Universe is accelerat- 
ing (Perlmutter et al. 1998, 1999; Schmidt et al. 1998; Riess 
et al. 1999) has led to a consensus that the cosmological 
model best fitting current data is a spatially-flat cold dark 
matter Universe with a matter density around 0.3 of the 
critical-density (Peebles 1984; Turner, Steigman & Krauss 
1984; Efstathiou, Sutherland & Maddox 1990). This model, 
known as ACDM, can boast an impressive range of obser- 
vational successes, with its main drawback being theoretical 
objection both to the magnitude and the required recent 
dominance of the cosmological constant term. 

It is often stated that while the supernova results are 
powerful in themselves, it is unlikely that they would have 
been widely accepted had there not been considerable other 
evidence pointing towards this favoured cosmology (Krauss 
& Turner 1995; Ostriker & Steinhardt 1995). Amongst that, 
one might mention the shape of the galaxy correlation func- 
tion, the combination of the cluster baryon fraction with 
standard Big Bang Nucleosynthesis (BBN), and the flat 
geometry inferred from the cosmic microwave background 
combined with the low matter density implied by direct ob- 
servations. 

In this paper, we examine the extent to which these 
additional arguments might be undermined in an alterna- 
tive cosmological model, which features an asymmetry in 
the lepton sector leading to a higher than usual abundance 



of neutrinos. It was claimed recently (Adams & Sarkar 1998, 
Lesgourgues & Peloso 2000) that these models offer one of 
the best remaining prospects for salvaging the idea of a 
critical-density Universe [the other main option being the 
Broken Scale Invariance models (Barriga et al. 2000)], and 
although the likelihood of doing so is small it is judicious to 
be aware of the possibility in order to balance its drawbacks 
with those of the cosmological constant model. 

The lepton asymmetry model relies on primordial pro- 
cesses to create an imbalance between the numbers of neu- 
trinos and antineutrinos in the Universe, which may reside 
in any of the three neutrino families. This would be the lep- 
tonic analogue of the (presently unknown) processes lead- 
ing to the baryon number of the Universe, though in this 
case interesting effects only arise for an asymmetry of order 
one, whereas the baryon-to-photon ratio is of order 10~^. 
There are many particle physics motivated scenarios for gen- 
erating such a large lepton asymmetry (e.g. Foot, Thom- 
son & Volkas 1996; Casas, Cheng & Gelmini 1999; March- 
Russell, Murayama & Riotto 1999; McDonald 2000; Dolgov 
et al. 2000; Kirilova & Chizhov 2000; Di Bari & Foot 2001). 
The lepton asymmetry leads to two important physical ef- 
fects. The first is that it modifies standard nucleosynthesis 
calculations, since the neutrino asymmetry alters the initial 
balance of protons and neutrons, and it has been known 
for some time that matching the element abundances in the 
presence of a strong lepton asymmetry can require a higher 
baryon fraction than standard nucleosynthesis (see for in- 
stance Kang & Steigman 1992; Esposito et al. 2000, 2001; 
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Kneller et al. 2001). The second is that it increases the ra- 
diation density in the Universe, by boosting the neutrino 
density beyond its usual value of 0.68 times the photon den- 
sity. 

At face value, these have highly desirable implications 
in cosmology for those favouring critical density on grounds 
of elegance, as was first pointed out by Adams & Sarkar 
(1998). First, if the preferred baryon density from nucleosyn- 
thesis could be significantly increased through the leptonic 
asymmetry, the cluster baryon fraction would then become a 
strong argument for critical density rather than against. Fur- 
ther, the extra radiation leads to a delay in matter-radiation 
equality, which shifts the characteristic bend in the matter 
power spectrum to larger scales mimicking the effect of the 
reduced matter density in the ACDM model. Finally, it was 
stressed (Lesgourgues & Peloso 2000; Esposito et al. 2001; 
Kneller et al. 2001) that a leptonic asymmetry could help 
in explaining last year's microwave anisotropy results from 
Boomerang (de Bernardis et al. 2000) and Maxima (Hanany 
et al. 2000), in which the observed weakness of the second 
acoustic peak favoured a high baryon fraction. Indeed, our 
initial studies for this present work indicated that critical- 
density models with leptonic asymmetry could fit not only 
these data, but also up-to-date constraints from large-scale- 
structure, the cluster baryon fraction, and primordial ele- 
ment abundances. So, although unable to explain the super- 
novae data, this model could undermine much of the other 
evidence supporting the ACDM model. 

This picture seems to be less promising after the 
publication of new microwave anisotropy results by DASI 
(Halverson et al. 2001; Pryke et al. 2001) and updated data 
analysis by Boomerang (Netterfield et al. 2001) and Maxima 
(Lee et al. 2001), which are incorporated into the results re- 
ported here. The new results contain no direct independent 
evidence for a cosmological constant, but within the frame- 
work of ACDM models exhibit excellent agreement with the 
baryon fraction from standard nucleosynthesis, while an ex- 
cess of baryons is a key ingredient for the success of critical- 
density models. However, a high baryon fraction may yet be 
allowed in the presence of a leptonic asymmetry, and the 
true test of the idea lies in detailed comparison with ob- 
servations, which is our purpose in this paper. We will see 
that the situation is not at all promising for the simplest 
case of massless neutrinos, which have trouble fitting many 
types of observation. However, there is now considerable ex- 
perimental evidence that neutrinos actually possess a small 
mass, and the effects of this need to be included. The neu- 
trino mass provides an additional modification to the matter 
power spectrum through neutrino free-streaming (as in the 
mixed dark matter scenario), and we find that this enables 
excellent fits to many observational data to be obtained. Un- 
fortunately, due to the latest Boomerang results, the model 
fails to explain the baryon fraction in clusters, as well as the 
present acceleration. 



and the mass of the neutrinos. In principle the masses at 
least are not extra parameters as compared to the standard 
cosmology, in that there is now substantial evidence that 
neutrinos do have mass; however in the presence of a lep- 
ton asymmetry the neutrinos may have a more significant 
impact on predictions for a given mass as the asymmetry 
increases their number density. In compensation for adding 
these extra parameters, we remove the cosmological con- 
stant. 

Provided that neutrinos reached thermal equilibrium 
before decoupling, the leptonic asymmetry for each flavour 
species can be conveniently parametrized by the ratio of 
chemical potential over temperature, ^j,. = fiv^/T^^ (with 
i £ {e, ^, r}). Neutrinos with a chemical potential are called 
degenerate neutrinos, because the asymmetry enhances their 
total density. When the neutrinos are in the relativistic 
regime, this effect is strictly equivalent to a change in the ef- 
fective number of standard neutrinos, in excess of the usual 
value of 3, of 

AiVeff = ^ [30(5.7^) V7 + 15(?.,/7r)V7] . (1) 

i 

All the generation mechanisms proposed so far predict 
different ^iz/s for each species, at least in absence of fine- 
tuning. This is a crucial point because the density of Ue and 
+ have opposite effects on the neutron-to-proton ratio 
at freeze-out during BBN, and on the production of light 
elements. More precisely, nucleosynthesis in the presence of 
a lepton asymmetry (known as degenerate Big Bang Nucle- 
osynthesis) requires three ingredients in order to be compat- 
ible with the observed abundances of deuterium, helium-4 
and lithium-7 : (i) an increas^ in i/e density (^^^ > 0); (ii) 
an increase in the baryon density; (iii) an increase in the to- 
tal density of radiation (and expansion rate of the Universe), 
bigger than the one resulting from (i), and parametrized by 
an effective number of standard neutrino species A'cff > 3. 
There are many possibilities for enhancing the radiation 
density: the mu and/or tau neutrino may have a leptonic 
asymmetry bigger than that of the electronic neutrino,^ or 
may become slightly non-relativistic during BBN (Hansen & 
Villante 2000), and apart from the three flavour neutrinos, 
many scenarios predict that extra relativistic degrees of free- 
dom could be present during nucleosynthesis (for instance, 
axions). Hannestad (2001) has recently studied limits on the 
number of neutrino species from the latest data (including 
models with a cosmological constant). 

Standard BBN, which corresponds to £,vi = 0, pre- 
dicts a baryon fraction given by Q.^ah^ = 0.019 ± 0.002, 
and an effective neutrino number close to 3. In the follow- 
ing analysis we will focus on a baryon fraction in the range 
0.015 < fib/i^ < 0.035. In this case, according to the most re- 
cent studies of degenerate BBN (Esposito et al. 2000, 2001; 
Kneller et al. 2001), the asymmetry parameter should 



2 THE LEPTON ASYMMETRY MODELS 

The lepton asymmetry models are in most respects the same 
£18 conventional structure formation models, in particular re- 
lying on the presence of cold dark matter, but add new pa- 
rameters describing the magnitude of the lepton asymmetry 



* Actually, there is also a small allowed region in parameter space 
with g,/^ < 0, reduced baryon density and N^g < 3, but this is 
irrelevant in the present framework. 

t i.e., 1^1,^ + I > i^y^. However, successful implementations 
require only a factor of order five or so between |^^^ -|- [ and 
^i/^ , which seems to be compatible with most mechanisms of large 
leptonic asymmetry generation. 
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be in the range < < 0.5, while the required effective 
neutrino number could be as big as 15 or even 20. Using 
equation we see that the contribution of ^i^^ to AA'^eff 
is negligible; therefore, when studying the spectrum of mi- 
crowave anisotropies and large-scale structure, we can forget 
completely about ^i,^, and consider the constraint from de- 
generate BBN to lie simply in the (Q.hh^ , Nes) plane. 

The calculations of matter and radiation power spectra 
were carried out using a modified version of the CMBFAST 
code (Seljak & Zaldarriaga 1996), as described in Lesgour- 
gues & Pastor (1999). 



3 OBSERVATIONAL CONSTRAINTS 

For each model, we define a statistic including the fol- 
lowing terms: first, 19 data points from the new analysis of 
Boomerang (Netterfield et al. 2001), 13 from the new analy- 
sis of Maxima (Lee et al. 2001), and 9 from DASI (Halverson 
et al. 2001); second, 22 data points from the PSCz redshift 
survey (Hamilton & Tegmark 2000); and finally, a constraint 
on the matter spectrum normalization as from the number 
density of galaxy clusters. For the last, we adopt the rather 
conservative constraint as — 0.56 ± 0.056 (l-f) from Viana 
& Liddle (1999). We also reran the analysis using the tighter 
limit as ~ 0.495±0.034 (l-a) recently obtained by Pierpaoli, 
Scott & White (2000), but this made no qualitative differ- 
ence to our conclusions and so we do not report those results 
here. 

For Boomerang and Maxima, we treat each data point 
as uncorrelated, with approximately gaussian window func- 
tions. We take into account the fully correlated calibration 
error and the multipole-dependent beam plus pointing error, 
minimizing over the corresponding parameters. The likeli- 
hood function can be written as C — e ^ ' provided that 
the data points are almost gaussian distributed. For Maxima 
we use the results for AT (Table 1 in Lee et al. 2001) which 
have almost symmetric errors (except for the last two points, 
which doesn't matter because these points provide mainly 
upper limits). The beam plus pointing uncertainty (calcu- 
lated from the same Table 1) turns out also to be symmetric 
for AT;, so we can define a for Maxima as 

^2 ^ {ATl^'° -{l + ca, + b a,j)AT[""^f ^ ^ ^ 
^—^ a? ' 

with a l-a calibration uncertainty ac = 0.04, and a l-a beam 
plus pointing uncertainty which is well fitted by the function 
Cb,i = 10"®/^ ''. For each model, we minimize the x^ over b 
and c. For Boomerang, we use a similar expression. However, 
in Netterfield et al. (2001), symmetric error bars are given 
for Di = {ATi f = 1{1-\-1)Ci/2tt. Therefore, we define the x^ 
directly on this quantity, with a l-a calibration uncertainty 
CTc = 0.20. The beam errors (read from figure 2 in Netterfield 
et al. 2001) are symmetric for AT;, with ab,; = 0.215x10"''/^ 
at l-a. For simplicity, we assume a gaussian beam error for 
Di, with twice the uncertainty. 

For the DASI data, Pryke et al. (2001) indicate that 
the use of the exact window functions, and of a transforma- 
tion that gives exactly gaussian errors (Bond, Jaffe & Knox 
2000), has only a modest impact on parameter extraction. 
On the other hand, the points cannot be treated as uncorre- 
lated. Accordingly, we define the following covariance matrix 



M,, = AA V^j AD, + s^D^D, , (3) 

using the data points Di ± AD; and the correlation matrix 
Vij from Halverson et al. (2001, Tables I and II). The fully- 
correlated uncertainty s equals 0.08. The x^ is then defined 

as 

2 X ^ / T-ithco T-,obs \ 1 ,f — 1 / 7-.thco 7-.obs \ ( t\ 

X =2^ (A -A )M^, (A -Dj ). (4) 

We compute the total values on a grid in parameter 
space, and perform a multidimensional cubic spline inter- 
polation in order to find the minimum and the confidence 
limits on each parameter.^ The precision and efficiency of 
this method depends crucially on the choice of a particu- 
lar parameter basis. Our seven cosmological parameters are 
the overall normalization (adjusted automatically to match 
the COBE observations (Bennett et al. 1996; Bunn & White 
1997) by the CMBFAST code that we use), the baryon frac- 
tion Q.hh^ , the scalar tilt tIs, the optical depth to reion- 
ization r, the mass of the degenerate neutrino rriv, the ef- 
fective neutrino number N^s (we recall that in the case of 
massive neutrinos, this number is defined at nucleosynthe- 
sis, not today), and a final parameter measuring the dark 
matter density. This last parameter could be taken as Q,A^h? 
(where f2dm = ^cdm+^v)', however, this choice would lead to 
an exceedingly large computing time because there is a de- 
generacy between Q.i^h? and the neutrino parameters (m^, 
Neff). In other words, for a given cosmological model and set 
of observations, the likelihood regions are elongated along 
a direction that can be found only empirically, and the 
varies very slowly when the function phi^dmh^ , Ncs,m,^) as- 
sociated with the degeneracy is almost constant. The best 
time-saving strategy is to directly use ph as the last free 
cosmological parameter. The preferred value of h (and of 
any other combination of the cosmological parameters) can 
then be recovered a posteriori. In most cases studied here- 
after, we find that ph = O.dmh'^ N~g {3.5 -\- m„)~^ is a fairly 
good parametrization of the degeneracy in the vicinity of 
the minimum (with expressed in electron-volts). 

In addition to the six free cosmological parameters, our 
model includes a free PSCz bias. The number of degrees of 
freedom is therefore (19 -I- 13 4- 9 -f- 22 -f 1) - (6 -I- 1) = 57. 
Actually, the constraints from PSCz on the largest scales are 
so loose that this number is somewhat overestimated. 

3.1 The massless neutrino case 

The results for massless neutrinos are summarized by the in- 
dividual parameter probability distributions, plotted in fig- 
ure |l|. The values for the best-fitting model and the 95 per 
cent confidence level are given in Table |l] 

t To marginalize over unwanted parameters, we maximize the 
likelihood function instead of integrating over these parameters 
(the two techniques would be strictly equivalent only for a multi- 
variate gaussian likelihood). Therefore, our confidence limit com- 
putation scheme is less rigorous than in current state-of-the- 
art analyses (Lange et al. 2000; Jafli'e et al. 2001; Netterfield et 
al. 2001), but because of its simplicity it is widely used by many 
other authors, and gives a fairly good hint of the true error bars 
[see for instance the discussion in Tegmark & Zaldarriaga (2000) 
and Tegmark, Zaldarriaga & Hamilton (2001)]. 
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Figure 1. The probability distribution for each cosmological 
parameter, in the case of one massless degenerate neutrino. The 
thick solid lines show the result including all the data. Parameter 
values are allowed at the 95 per cent confidence level when the 
probability exceeds the horizontal line. The thin curves show the 
constraints obtained by combining just one CMB experiment with 
the other non-CMB data: solid is Boomerang, dashed is Maxima 
and dot-dashed is DASI. Although Maxima favours a significantly 
higher baryon fraction, and Boomerang a lower scalar tilt, the 
three data sets are found to be perfectly compatible. 

The results for this model are quite disappointing. The 
best-fitting model has a of 61, which given the number 
of degrees of freedom looks quite satisfactory. However the 
properties of the best-fitting model are quite undesirable. 
The preferred values f2b/i^ = 0.018 and A'ofi = 3 match the 
standard BBN prediction, but this is not a good thing in the 
present context; since we don't obtain a high baryon density, 
we cannot explain the cluster baryon fraction (our predicted 
value for 57b^^'^ and the observed lower bound (Ettori & 
Fabian 1999) have no overlap at the 2-a level). A further 
problematic aspect of this model is its low preferred value of 
h: this gives an impressively large age, but is in considerable 
discrepancy with direct h measurements. Finally, we must 
recall that we have made no attempt to obtain a presently 
accelerating universe. 

We are therefore forced to conclude that the massless 
neutrino case has too many failings against observations to 
be considered a viable model. 

3.2 The massive neutrino case 

We now suppose that the neutrino family with the lep- 
tonic asymmetry^ has a mass m^. Now the degenerate neu- 

§ In the massless case our results were model-independent, since 
we did not privilege a particular scenario for the origin of the 
extra relativistic degrees of freedom. When taking into account a 
neutrino mass niy, we could distinguish various cases: first, the 
large chemical potential responsible for Af(,ff > 3 during BBN 
can belong to the massive neutrino family; alternatively, it can 
be shared between one species with negligible mass and one with 
arbitrary mass mu\ finally, the extra radiation density could be 
attributed to particles other than flavour neutrinos. For brevity, 
we only discuss the simplest case of a single massive degenerate 
neutrino family (i^^ or i^^)- Most other situations would give com- 
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Table 1. The preferred value and 95 per cent confidence limits 
for each cosmological parameter, in the case of one massless and 
one massive degenerate neutrino family. The upper part refers to 
the parameter basis used in the interpolation, with ph defined as 
in the text. The lower part refers to useful combinations of these 
parameters, including the leptonic asymmetry parameter ^u, the 
age of the Universe to and the quantity Q^,h^-^' which can be 
compared with the baryon fraction in galaxy clusters (the range 
given by Ettori & Fabian (1999) is 0.060 ±0.025 at 2-a confidence 
level) . 

trino can make up a significant fraction of the dark mat- 
ter, as in the mixed dark matter scenario, and its free- 
streaming while relativistic suppresses short-scale matter 
perturbations. With this additional free parameter, the min- 
imum of shows a large degeneracy along ph ■ unreasonably 
large values of h are allowed, with a huge efi'ective neu- 
trino number maintaining the first acoustic peak and the 
power spectrum with the right shape and amplitude. This 
parameter region is uninteresting and should be removed. 
Indeed, h > 0.58 corresponds to a Universe younger than 
to = 11 Gyr, which is almost completely excluded. So, we 
must add to the ^ "weak age prior" to > 11 Gyr It is 
important to note that this prior almost does not affect the 
goodness-of-fit of the model, since the best-fitting model has 
to close to 11 Gyr anyway. 

The results for the massive case are also given in Ta- 
ble 0, and the probability for each individual parameter is 
shown in figure |2[ The best-fitting model now has an impres- 
sively low of 43 and some remarkable features. A large 
effective neutrino number between 7 and 15 is preferred, 
producing a high first acoustic peak as in ACDM models. 
This large lepton asymmetry is compatible with BBN up 
to A^cff — 9, as can be seen in figure ^. When it is com- 
bined with a mass close to 1 eV, it gives the right shape and 
amplitude for the matter power spectrum. A neutrino mass 
smaller than 0.6 eV is excluded at more than 95 per cent 

parable results, but with a higher preferred value of the mass, 
since the neutrino free-streaming effect is enhanced by the lep- 
tonic asymmetry (Lesgourgues & Pastor 1999). 
^ Technically, this is done by multiplying the likelihood function 
by a gaussian cut-off, if and only if to < 11 Gyr. The variance is 
chosen so that at to = 10 Gyr the cut-off factor equals 1/2. 
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Figure 2. The probability distribution for eacli cosmological 
parameter, witli one family of massive degenerate neutrinos. As 
in figure hi the thick solid line shows the results from the com- 
plete data set, while the thin curves combine the non-CMB 
data with each of the three CMB experiments individually. Max- 
ima and DASI allow a significantly higher baryon fraction than 
Boomerang, though the three data sets are compatible. 



confidence; this result is in nice agreement with the osciUa- 
tions reported at the Los Alamos Liquid Scintillation Neu- 
trino Detector (Athanassopoulos et al. 1998), which support 
evidence for a neutrino mass > (0.1 — 1) eV^. 

Unfortunately, this positive picture is darkened by the 
predicted baryon density, which is as low as in the standard 
case: Qhh^ = 0.0201q qq4 (95 per cent confidence). So, within 
the range of viable parameters the lepton asymmetry model 
can no longer reach the high baryon fraction potentially al- 
lowed by the degenerate BBN model. Studying the curves for 
individual CMB experiments in figure ^, we see that this re- 
sult is driven primarily by the new Boomerang results, with 
DASI and Maxima both still allowing significantly higher 
values; note in particular that inclusion of neutrino mass 
allows DASI to go to higher baryon fractions than it can 
in the massless case. With all data taken into account, the 
model is now restricted to fib/i'^'^ < 0.033, more than 2-a 
away from the Ettori & Fabian (1999) cluster bound. 



4 SUMMARY 

We have performed a detailed comparison of critical-density 
models including leptonic asymmetry with the latest CMB 
and LSS data, in order to investigate their viability as al- 
ternatives to the ACDM model. Some sample power spectra 
are shown in figure ^. We have found that very good fits to 
those data are available, due to the combined effect of the 
large lepton asymmetry (which is compatible with primor- 
dial abundances), and of a neutrino mass ~ 1 eV (which is 
in nice agreement with LSND). This model cannot hope to 
explain the supernovae data, but has the prospect of under- 
jcha btha support for the ACDM paradigm which h-AH 
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data introduces a new problem for this model, which is that 
the baryon density is now constrained to be low enough that 
fits to the cluster baryon fraction are not possible, which is 
disappointing as the lepton asymmetry model had the po- 
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Figure 3. The thick lines show the allowed region in the {fl^Jh^, 
N^ff) parameter space, at 68 and 95 per cent confidence levels, for 
one family of massive degenerate neutrinos. On the lower axis we 
show the standard BBN prediction, and the thin curves show the 
region allowed at 95 per cent confidence level by degenerate BBN 
(Esposito et al. 2000, 2001). The regions overlap for 5 < N^s < 9. 



tential to allow higher baryon densities while remaining com- 
patible with nucleosynthesis. The main driving force to this 
conclusion is the new analysis of the Boomerang data (Lee 
et al. 2001); the other new CMB data still permit a higher 
baryon density in the presence of a massive degenerate neu- 
trino. Given the subtle effects of the neutrino degeneracy, 
and the various uncertainties in the new CMB data (cali- 
brations, tilts), this was not obvious by eye, and it is the 
main result of this paper. 

We stress that our study does not provide a model- 
independent bound on the leptonic asymmetry in the Uni- 
verse, since it could in principle coexist with a cosmological 
constant. However, the lepton asymmetry is better moti- 
vated in the critical-density case, with it being used to re- 
move the need for A. What our study shows is that follow- 
ing the recent results, the critical-density lepton asymmetry 
model experiences new observational difficulties which make 
it a less attractive proposition as a simple and elegant alter- 
native to the ACDM cosmology. 
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Figure 4. The CMB anisotropy and matter power spectra for three models: the massless and massive degenerate neutrino best-fitting 
models (whose parameters are given in Table hi), and the preferred ACDM model given by Wang, Tegmark & Zaldarriaga (2001): 
n^h^ = 0.020, fl^dmh'^ = 0.012, Oa = 0.66, = 0.93. The corresponding (including CMB and LSS data) are equal to 61, 44 and 
56 (the last is not very good simply because Wang et al. did not include a ag constraint in their analysis). In the first three graphs, the 
Boomerang, DASI and Maxima data sets are shown with the appropriate beam and calibration errors (6, c) calculated for each case. The 
three matter power spectra are plotted together in the final plot, along with the PSCz points divided by the square of the bias factor 
fe = 1.4 which minimizes the for the massive neutrino model. 
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